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COMPUTATION OF TOPSIDE lONOGRAMS FROM N(h) PROFILES 


By Lawrence Colin and Kwok-Long Chan 
Ames Research Center 
and 

Jack G. K. Lee 
Informatics, Inc. 

SUMMARY 


Model profiles of plasma frequency and gyrofrequency are chosen to 
simulate typical conditions observed in the topside ionosphere by the Alouette 
I sounder (1000 km) and the Alouette II sounder near apogee (3000 km) . Exam- 
ples of group refractive index profiles for these models are presented. By 
means of parabolic-in- log X lamination techniques applied to the plasma fre- 
quency profiles, ’’cumulative delays” are computed to illustrate the relative 
delays contributed by each lamination to the formation of topside ionograms . 

"Complete ionograms (ordinary, extraordinary, and Z traces) are presented for 
all the models. Implications of the results to the inverse experimental prob- 
lem of scaling ionogram data and computing electron concentration profiles are 
discussed. 


INTRODUCTION 


Whereas research concerned with ionospheric sounding from the earth’s 
surface is more than 25 years old, sounding of the topside ionosphere from 
satellite-borne radars has been conducted for only six years. The satellites 
Alouette I (1962 to present, circular orbit at about 1000 km) and Alouette II 
(1965 to present, elliptical orbit with apogee at 3000 km, perigee at 500 km), 
containing HF swept- frequency sounders, have provided the bulk of available 
data on electron concentration, N(h), above the F 2 layer peak. 

Because of the complexities of topside ionograms, sophisticated N(h) 
reduction procedures and programs have been developed. The standard reduction 
procedure involves computing the N(h) profile from a set of virtual depth- 
frequency pairs, h*(f), scaLed from the extraordinary trace appearing on a 
filmed ionogram. This technique is adequate for the accurate reduction of the 
great bulk of Alouette I and low-altitude Alouette II ionograms (ref. 1) . On 
the other hand, there remain a great number of topside ionograms, particularly 
those collected at the higher orbital altitudes of Alouette II, for which 
the standard reduction procedure yields inaccurate, often unrealistic, elec- 
tron concentration profiles. The difficulties are caused mainly by scaling 
uncertainties associated with large virtual depths, h’(f), steep gradients. 



dh’/df, sharp "cusps,” and the occurrence of reflection traces 

arising from oblique propagation. Clearly, new procedures involving optimum 
use of independent and/or redundant information available in the ionogram data 
must be brought to bear if reliable high-altitude N(h) profiles are to be 
computed. For example, the use of the redundant information available in all 
reflection traces may well enhance the accuracy of computed profiles. 

The present study deals with the computation of ionograms from model 
N(h) profiles. This direct approach clearly demonstrates the transformation 
performed by the ionosphere on the sounder waves and the one-to-one corre- 
spondence of true and virtual reflection heights (ref. 2) . The results supply 
insight to the redundancy problem and suggest some feasible schemes for opti- 
mum implementation of N(h) reduction procedures using extraordinary, 
ordinary, and Z trace data. 


REFRACTIVE INDICES IN THE TOPSIDE IONOSPHERE 


The transformation is described by 


h’(f) = p* dh (1) 

which relates the virtual depth, h*, of reflection of a particular propagation 
mode at a sounder frequency, f, to the true height of reflection, h^, for a 
given profile of group refractive index, y * , between the altitude of the sat- 
ellite, hg, and h^. Appendix A details the appropriate group refractive index 
and phase refractive index, y, equations expressed in terms of the familiar 
magnetoionic parameters X, Y, and 6. Before proceeding with the major pur- 
pose, it is instructive to examine carefully the variations of y and y^ with 
these parameters. 

The variation of y^ and y^ with X is illustrated in figure 1 for a 
value of Y < 1 and in figure 2 for a value of Y > 1. Curves are shown for 
several values of 0. A combination of linear and reciprocal scales is used 
to include both zero and infinity on abscissa and ordinate. For Y < 1 
(fig. 1) there is one branch for the ordinary mode (0 ^ X < 1) and two 
branches for the extraordinary mode (0 < X < 1 - Y and 

(1 - Y2)/(1 - Y2 cos^ 6) < X < 1 + Y) . The latter branch of the extraordinary 
mode is referred to as the Z mode. For Y > 1 (fig. 2) the ordinary mode 
has two branches (0 < X < 1 and (1 - Y^)/(l - Y^ cos 9) < X; the latter branch 
is called the whistler mode and exists only when Y^ cos^ 0 > 1) while the 
extraordinary mode has only one branch, the Z mode branch (0 ^ X < 1 + Y) . 
The ordinary, extraordinary, and Z modes are consistently observed as dis- 
tinct reflection traces in topside ionograms . Nonzero values of X at the 
satellite sounder antennas allow the observation of the Z mode for Y < 1. 

In bottomside ionograms the Z trace is observed for Y < 1 only when wave 
coupling exists between the ordinary and extraordinary modes (ref. 3) . 
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Zeros in occurring at X = 1 (ordinary mode) , X = 1 - Y (extra- 
ordinary mode) , and X = 1 + Y (Z mode) correspond to wave reflection condi- 
tions for the respective modes . An infinity in occurs at 

X = (1 - Y^)/(l - y2 cos^ 0) (Z mode for Y < 1, and ordinary mode for 

Y > 1) . Clearly, both zeros and infinities in y^ correspond to infinities 

in y*. However, the infinity in y’ associated with 

X = (1 - Y^)/(l - Y^ cos^ 9) is not a wave reflection condition as are the 
other infinities. 

Graphs of y^ and y* similar to figures 1 and 2 have been published in 
references 3-6. In topside N(h) reduction, it is necessary to take into 
account the Y variation with height. Thus it is useful to consider the varia- 
tion of y^ and y’ with X and Y for a fixed value of 6, since 0 is 
nearly independent of height for vertical propagation. Figures 3 and 4 illus- 
trate the variations of y^ and y’ versus X for 0=5° and 55°, respectively, 
and for several values of Y for each propagation mode. These particular 
cases of 0 are used in the analyses reported below. NASA TM X-1553 (ref. 7) 
contains many graphs of the types illustrated in figures 1 to 4. 


TOPSIDE MODELS 


Common types of ionograms observable from the Alouette sounders have been 
computed from equation (1) using several realistic model profiles of plasma 
frequency, fj^(h), and gyrofrequency , f|^(h), shown in figure 5. High and low 
f^(h) models and high and low fh(h) models extend from h = 3000 km to 400 km. 
The two fn(h) models have identical shapes corresponding to a theoretical 
plasma distribution based on vertical diffusive equilibrium such that the 
electron concentrations at any altitude differ by a factor of 5 between models 
(ref. 1). The high and low fh(h) models are dipolar (inverse cube) with mag- 
nitudes at 3000 km typical of those observed at high latitudes and low mid- 
latitudes, respectively. The propagation angles, 0, are independent of h 
with magnitudes 5° and 55° to simulate vertical propagation at high latitudes 
and low midlatitudes. A total of eight simulated Alouette cases, summarized 
in table I, were constructed from these models. 


^’CUMULATIVE DELAY” CURVES AND IONOGRAMS 


For each case in table I, ionogram traces for each propagation mode may 
be calculated by integrating equation (1) using the group refractive index 
equations given in appendix A. The integration procedures are described in 
appendix B. The fn (h) and fh(h) profiles are subdivided into equal thickness 
contiguous laminations (100 km for cases 1 through 4, 40 km for cases 5 
through 8). Within each lamination the height is assumed parabolic-in-log X 
(see appendix B) and is assumed to vary in an inverse-cube manner. Using 

the reflection criteria described above, that is, ordinary mode (X=l) , extra- 
ordinary mode (X=l-Y) , and Z mode (X=l+Y) , those frequencies reflected at 
the lower boundary of each lamination are calculated. These sets of 
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frequencies are used to compute the ionogram traces. Six examples are 
illustrated in figures 6 through 11. Part (a) of each figure shows the group 
refractive index profiles, y'(h), appropriate to each lower boundary reflec- 
tion frequency. The corresponding ionogram trace is given in part (b) where a 
smooth curve is drawn through the computed virtual delays for each reflection 
frequency. Also shown in part (b) are "cumulative delay" curves which illus- 
trate the relative delays contributed by each lamination as the sounder wave 
"propagates" to the reflection level. 

Figure 6 represents computation of the extraordinary mode for case 1. At 
each frequency, y’Ch) increases monotonically and smoothly from a finite value 
at the satellite to infinity at the altitude of reflection corresponding to 
the condition X = 1 - Y, Integration of the y’(h) curves produces a moder- 
ate "cusp" in the extraordinary trace at frequencies slightly greater than the 
satellite cutoff frequency, fxs • The "cumulative delay" curves quite closely 
parallel the trace, and the delay in the lowest lamination, for each frequency, 
is significantly greater than the delays in each of the preceding laminations. 
These variations are typical of extraordinary mode variations for all eight 
cases. Subtle differences do occur between the results of the high and low 
fh(h) models because the extraordinary mode group refractive index is sensi- 
tive to Y, but not to 0, as seen in figures 3 and 4. 

Figures 7 and 8 represent two examples for the ordinary mode (cases 1 
and 2) in which the y*(h) profiles are quite similar to each other, that is, 
y’(h) varies smoothly from finite values at the satellite to infinity at those 
altitudes where X = 1. However, the computed ionograms are quite dissimilar. 
Figure 7 illustrates the formation of a sharp "cusp" caused by an enormous 
delay encountered in the lamination of wave reflection which is typical for 
high-latitude ionograms where 0 is small. Figure 8 illustrates the forma- 
tion of a more moderate cusp. Although the delay in the last lamination is 
still largest relative to each of the preceding laminations, the absolute 
delay is considerably smaller than that in figure 7. The contrasts seen in 
figures 7 and 8 demonstrate that the ordinary mode group refractive index is 
a sensitive function of 0, but not of Y, as shown in figures 1 and 2. 

The y*(h) profiles and ionograms for the Z mode for cases 2, 4, and 8 
are shown in figures 9, 10, and 11. Figure 9 shows that for low frequencies 
the variations of y*(h) are similar to those in figures 6, 7, and 8. How- 
ever, at higher frequencies (f > 0.228 Mc/s) a nonmonotonic behavior occurs 
when y’ (h) decreases with depth to a minimum value before increasing to 
infinity at the reflection level where X = 1 + Y. Finally, at a frequency, 
fzoo? corresponding to X = (1 - Y^)/(l - Y^ cos^ 0), y' becomes infinite at 
the satellite and the Z mode ceases to propagate from the sounder antennas. 
This behavior is most frequently observed on topside Z traces. However, it 
is not unusual to encounter the situations shown in figures 10 and 11. Here, 
for low and medium frequencies, the y*(h) profile behavior is quite similar 
to that described for figure 9. However, at frequencies (f > 0.290 Mc/s in 
fig. 10 and f ^ 0.726 Mc/s in fig. 11) close to fzoo> "^he y’(h) profiles 
take the peculiar variations shown, that is, y’(h) exhibits two minimums and a 
maximum before increasing to infinity at the reflection level. This behavior 
is usually encountered for conditions of high fh (h) and low fn(h) where Y 
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is greater than 1 and X varies from less than to greater than 1. The well- 
defined maximums in figure 10 are due to the maximum of near X = 1 

observed in figure 4(b) for values of Y near 2. The large values of p’(h) 
immediately below the satellite height for f = 0.959 Mc/s in figure 11 are 
due to the large values of observed in figure 4(b) in the region near and 

below X = 1 for values of Y only slightly greater than 1. These peculiar- 
ities in the y'(h) profiles give rise to the interesting Z traces seen in 
figures 9, 10, and 11. To produce the detailed variations seen in these fig- 
ures, lamination thicknesses of 20 km had to be used to provide adequate fre- 
quency resolution. However, only the 100 km (fig. 9(a)) and 40 km (fig. 10(a)) 
"cumulative delay" curves are shown. The Z trace with a smooth, monotonic 
variation of virtual depth with frequency, shown in figure 9(b), is most fre- 
quently observed in topside ionograms. On the other hand, figure 10(b) shows 
the development of a "plateau" in the Z trace. At low frequencies the vari- 
ation is smooth, as in figure 9(b), with the largest delay in the lamination 
of wave reflection. At f near 0.29 Mc/s, where the y’(h) maximum first 
occurs (shown as a dotted line) , a precipitous increase in virtual depth is 
observed. For 0.29 Mc/s < f < f^^, the virtual depth remains relatively con- 
stant. In this region the largest relative delay is provided by that lamina- 
tion exhibiting the y’(h) maximum. As f approaches 'the normal rapid 

increase in virtual depth is again observed with the largest delay occurring 
in the lamination immediately below the satellite. These features, although 
not quite as drastic, are seen in the Z trace in figure 11 simulated for an 
Alouette I ionogram. 

These peculiar Z traces occur at rather low frequencies where the 
Alouette sounder response is weak and at high latitudes where ionograms are 
normally poor. Also, the delays become quite large and some of them would 
only be observable on a "second-time-around" basis. Thus actual observation 
of these shapes is relatively infrequent. With improved sounder design for 
future satellites, however, observation should be possible and its use for 
N(h) reduction become feasible. 

The entire set of computed ionograms (ordinary, extraordinary, and Z 
traces) for all eight cases are shown in figures 12 through 19. The low- 
altitude ionograms are plotted on a linear frequency scale simulating Alouette 
I ionograms. The high-altitude ionograms are plotted on two different linear 
frequency scales for below and above 2 Mc/s, simulating Alouette II ionograms. 


DISCUSSION 


Figures 12 through 19 illustrate most types of ionogram traces due to 
vertical propagation which are observed in Alouette I and Alouette II sounder 
data. At h’ = 0 the traces always occur such that f^s > ^os ^ ^zs s^^id 
(dh^/df) = (dho/df) = dhz/df = «>. Although the cusps are often observed on 
ionograms from both Alouette I and Alouette II, they appear only in the simu- 
lated high- altitude ionograms shown in this study. It can be shown that the 
plasma scale height, defined by -N(dN/dh)”^, must have a height gradient near 
the satellite that exceeds a certain critical value for the cusp to appear 
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(J. Shmoys, private coiranuni cation) . Because the model profiles are not 
defined below 400 km, F 2 layer penetration is not simulated on the computed 
ionograms . The ionograms shown were computed for profiles having lamination 
thicknesses of 100 km and 40 km for Alouette II and Alouette I, respectively, 
corresponding to a total of 27 and 16 ’^scaled” points (or 26 and 15 
laminations) . 

It is difficult, in practice, to prescribe ionogram scaling criteria for 
optimum N(h) computation. It is, in principle, desirable to scale the iono- 
gram trace so that the computed N(h) profile has laminations of equal 
thickness, say, of about 100 km for Alouette II and of about 50 km for Alou- 
ette I data. This would require a total of about 25 points for the former and 
15 points for the latter ionograms . Computation times on large computers are 
then typically a few seconds, which is reasonable for the vast number of top- 
side ionograms available for study. The **equal lamination thickness” scaling 
criterion is impossible to apply strictly since, a priori, the N(h) profile 
is unknown. However, the models studied do illustrate that, for this criter- 
ion, the scaling procedures are quite different for the cusp and no-cusp cases. 
For no- cusp, the scaled points are separated by approximately equal intervals 
along the traces. For the cusp, the scaled points are thinned out in spacing 
where the trace curvature is small and packed in spacing where the curvature 
is large - the greater the curvature, the closer the packing. 

Although the simulated high-altitude Alouette II ionograms (figs. 12-15) 
and Alouette I ionograms (figs. 16-19) have many similarities, they have 
subtle differences that are observable in both simulated and actual ionograms. 
The Alouette II sounder frequency sweep rate below 2 Mc/s is about six times 
slower than the Alouette I sweep rate. Thus, greater frequency resolution is 
available in the regions where the large slopes occur. Although in principle 
the virtual depth resolution is identical for the two sounders because of 
identical pulse widths, in practice the resolution is poorer on Alouette II 
ionograms since the virtual depth range coordinate, of about 4500 km, occupies 
the same film dimension as does the 1500 km virtual depth range for Alouette I. 
Therefore, the available h’(f) resolutions are about the same for Alouette I 
and II ionograms. However, the plasma scale height increases with altitude 
in the topside ionosphere, and the Alouette II trace may have much larger 
virtual depths and slopes immediately above fxs^ ^ os > ^zs* available 

resolutions which are enough to yield accurate computed N(h) profiles for 
most of the Alouette I ionograms are not sufficient for many high-altitude 
Alouette II ionograms. Much more accurate and careful analysis is required. 

Extraordinary traces in topside ionograms are usually continuous and 
complete while the ordinary trace generally exhibits gaps between the plasma 
frequency (or the gyrofrequency when f^g ^ ^n) upper-hybrid frequency, 
due to antenna impedance mismatch, and the Z trace is usually weak, due to 
noise, interference and low sounder sensitivity, and/or excessive dispersion. 
Thus, the extraordinary trace is generally used, and for Alouette I is the 
only trace which has been used, for routine N(h) reduction. In actual iono- 
grams a single reflection trace for each propagation mode is normally observed 
and, in most cases, the traces are due to vertical propagation. However, 
traces due to oblique propagation may appear with or without the appearance of 
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a vertical propagation trace. For high-altitude Alouette II ionograms, 
because o£ the scaling accuracy difficulties and because of the frequent 
appearance of multiple traces, particularly for the extraordinary mode, it 
would be desirable to include the ordinaiy and Z trace data in the reduction 
procedures. Comparison of figures 12 through 19 shows that the general shape 
of the extraordinary trace remains quite stable from case to case, certainly 
more so than the ordinary and Z traces . As mentioned earlier, the extra- 
ordinary group refractive index is a more sensitive function to Y than 0, 
while the ordinary group refractive index is more sensitive to 0 than Y. The 
variation of Y (i.e., f]^) is quite small during the passage of the satellite 
over one telemetry recording station whereas 0 may change significantly. 

The sensitivity of the ordinary trace to 0 is further illustrated in fig- 
ure 20. Here the computed ordinary traces for case 1 are shown for the same 
fn and fh model profiles but with 0 allowed to vary from 5° to 30°. Hence, 
this difference in sensitivity of each mode to the magnitude and angle of the 
magnetic field vector indicates that the use of redundant data may improve the 
accuracy of computed N(h) profiles. 

An operationally feasible closed- loop reduction procedure that uses 
redundant data to improve computed N(h) profile accuracy and/or operator con- 
fidence is summarized as follows: (1) Scale the extraordinary trace; (2) com- 
pute an N(h) profile therefrom; (3) compute ordinary and Z traces from the 
computed N(h) profile; (4) compare the computed ordinary and Z traces with 
the visible portions of those traces on the ionogram; (5) if necessary, cor- 
rect the extraordinary trace scaled points and repeat this procedure until the 
comparisons in step (4) agree. Since the satellite moves in position while 
the traces are being formed this comparison becomes questionable should the 
N(h) distribution drastically change during this period. 

Another useful procedure is to (1) scale the extraordinary and Z traces 
with a least-squares fitting technique to deduce an accurate N(h) profile 
near the satellite; (2) from this partial profile compute the ordinary trace 
to bridge the gap where that trace is not observed; (3) continue the N(h) 
reduction using the ordinary trace, in conjunction with the extraordinary 
trace if necessary. The particular advantages of this approach are that 
(a) the ordinary trace is normally spread and dispersed much less than the 
extraordinary trace and (b) the ordinary trace more likely results from 
vertical propagation. 

Yet another approach would be to (1) scale the visible portion of the 
three traces; (2) compute a profile from each; and (3) combine the three 
individual profiles in a least-squares manner to compute a final N(h) profile. 

For those high-altitude ionograms in which the visible traces are all due 
to oblique propagation, the above procedures will not yield accurate profiles. 
Then, ray tracing analyses must be performed to determine the ionospheric 
distribution. 
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A data reduction system, called FILMCLIP, is currently operational at 
Ames Research Center and incorporates the features mentioned. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, May 6, 1968 
188-39-01-01-00-21 
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APPENDIX A 


REFRACTIVE INDEX EQUATIONS 


The group refractive index defined by 

y' = c/Vg 

is obtained from the phase refractive index 

y E c/vp 

by the process 

y’ = d(yf)/df 

where 

c velocity of light in free space 

Vg group velocity of the electromagnetic wave 

Vp phase velocity of the electromagnetic wave 

f sounder frequency 

The collisionless phase refractive index is given by the Appleton 
equation of the magnetoionic theory (ref. 3): 


- 1 1/2 


sin^ Q 
2(1 - X) 


I 


Y4 


4(1 


sin*^ 0 

w 


+ Y2 


COS'^ 


where 

X = N/12,400f2 

N electron concentration, electrons/cm^ 

£j^ electron plasma frequency, Mc/s 

Y fh/f =2.8 B/£ 

B magnetic induction o£ earth's field, gauss 

f}^ electron gyrofrequency , Mc/s 


(Al) 

(A2) 

(A3) 

Hart re e 
(A4) 
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angle between the wave normal and earth's magnetic field directions 
(0 - 90° - I, where I is the magnetic dip, for vertical propagation 
in a stratified ionosphere) 

It can be seen that two modes of propagation are possible in a magnetoionic 
medium, the ordinary mode (with the plus sign) and the extraordinary mode 
(with the minus sign). These are usually denoted with appropriate subscripts 
X and o. For clarity, the ordinary and extraordinary phase refractive indices 
may be written: 





(AS) 


Fx 



where 


sin^ d YR 

2(1 - X) 2(1 - X) 


(A6) 


(A7) 


Y^ sln^ 6 YR 

2(1 - X) 2(1 - X) 


(A8) 


R = [Y^ sin^ 


0+4 cos2 0(1 - 



(A9) 


Applying equation (A3) to equations (AS) through (A9) yields the ordinary and 
extraordinary group refractive indices: 


M-o 


XY(l - X)cos2 0 
S^(Y sin.2 0 + r) 


1 - 


Y sin^ 0 fl + 


R 


-S]} 


(AlO) 


X Hx I 


X(1 - Sx) 


2S| 


1 + 


Y sin^ 0 /l + X 


R 


1 - X 


(All) 


These equations simplify considerably for the special cases of transverse 
propagation (0 = 90°) and longitudinal propagation (0 = 0°) : 
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(a) Transverse propagation 


= (1 - X) 


1/2 




1 - 


X(1 - X) 

1 - X - Y2 


11/2 


I 1 


^"x 


1 + 


XY^ 


l^x L (1 
(b) Longitudinal propagation 


- X - Y^) -I 


= 1 - — ^ 
o V 1 + Y 


1/2 


^^x = (1 - 


1 - X 


1/2 


' _ 1 


1 - 


XY 


2(1 + Y)^J 


' _ 1 
" l^x 


1 + 


XY 


2(1 - Y)‘ 


(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 


Both Px 1‘x longitudinal case exhibit discontinuities at 

X = 1 when Y > 1, that is, where the extraordinary mode is generally known 
as the Z mode. It can be shown (ref. 4) that the following relationships 
result : 


^^x = 



for X < 1 

for X > 1 


(A20) 

(A21) 
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"i + ^ 1 


2(1 - Y)^J 


[l - ^ ' 

M*z 1 

L 2(1 + Y)^_ 


for X < 1 

for X > 1 


For 9^0° and Y > 1, exhibits a maximum at X = 1 and can be 

approximated by 




t 

z 


s 1 + 


Y2 


_1 

sin^ 6 


(A22) 

(A23) 

(A24) 
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APPENDIX B 


COMPUTATIONAL PROCEDURES 


Appendix A summarizes the equations for y’ as functions of X, Y, and 
0. For any frequency, f, profiles of X and Y are simply computed from the 

given fn(h) and fh(h) models, and 0 is assumed height invariant. The 

heights of reflection, h^, for the ordinary, extraordinary, and Z modes are 

given as the heights where X = 1, 1 - Y, and 1 + Y, respectively. Hence, all 

quantities on the right-hand side of equation (1) are known and the integra- 
tion for virtual depth, at that frequency, is easily calculated. In actual 
ionogram analysis where hp is the unknown, given h*, the situation is more 
complex. Since the parameter Y appears in the integrand and its integral 
upper limit, an iteration scheme is required to insure convergence to the 
correct reflection height (ref. 1) . 

The cumulative delay curves and ionogram traces were calculated for the 
fn(h) and f^Ch) models by numerically integrating equation (1) on a digital 
computer. The fn profiles were divided into 100 km and 40 km laminations 
for the high-altitude (3000 km) and low-altitude (1000 km) profiles. Equa- 
tion (1) was then solved for those frequencies for each propagation mode hav- 
ing reflection heights that occur at the lower boundary of each lamination. 
Within each lamination the parabolic-in- log X assumption was used, that is, 
in the (i + l)st lamination: 

h = hi + ai+i log ^ + bi+^log ^ (Bl) 

where hj^ and Xj^ are the values at the upper boundary of that lamination. 

The gyrofrequency was assumed to have an inverse- cube variation between its 
values at the upper and lower boundaries. Continuity in slope was assured at 
each boundary by the condition: 

^i+i = + 2h^ log ^ (B2) 


A parabolic- in- log X starting procedure was also used. Since no boundary 
condition is available at the sounder, however, the first two laminations are 
used to define the initial parabolic variation. 


Applying these lamination techniques, equation (1) becomes 




dh ; 


M 




N - 1 


(B3) 


and 
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^ n^±+l 

h'(f) = ) / (B4) 

i=i 

where is "cumulative delay," N - 1 is the lamination number, and h^ = hg . 

Since y' becomes infinite at the reflection heights (and at the sounder 
for fz^) , special care must be taken for accurate integration. Changes of 
variables are made in those lamina where y* becomes infinite such that the 
modified integrand is finite. These variable changes are given by: 

tQ = 1 - X for all laminations in the ordinary mode 



for all laminations in the extraordinary mode 


tf = 1 r for the first lamination in the Z mode, and 

^ (1 - cos2 


tf = 1 — for all but the first lamination in the Z mode 

^ 1 + Yr 

where is the value of Y at the height of reflection. Gaussian quadra- 

ture is used to evaluate the modified equation (B3) and (B4) integrals. It 
has been shown that, in general, three-point gaussian coefficients provide 
sufficient accuracy (ref. 1) . However, for small 0 higher order gaussian 
quadratures are required for the ordinary mode. In general, it suffices to 
use 


3-point gaussian for 0 > 40^ 
7-point gaussian for 10° < 0 < 40° 
16-point gaussian for 0 < 10° 
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TABLE I.- PLASMA AND GYROFREQUENCY MODELS 


ON 


Case 

hg, km 

Description 

0, deg 

f„s, Mc/s 

Ng , cm ^ 

fhs- 

Simulation of 

1 

3000 

high 

high 

fh 

5 

0.64 

5x103 

0.50 

Alouette II apogee, 

high latitude, day | 

1 

2 

3000 

low 

low 

fh 

55 

0.28 

103 

0.30 

1 

Alouette II apogee, | 
low latitude, night 

3 

! 

t 

3000 

high 

£^, low 

% 

55 

0.64 

5x103 

0.30 

Alouette II apogee, 
low latitude, day 

i 4 

1 

1 

1 

3000 

low 


fh 

5 

0.28 

103 

0.50 

Alouette II apogee, 
high latitude, night, 

5 

1000 

high 

high 

^h 

5 

1.49 

2.8x10^ 

1.03 

Alouette I, high 
latitude, day 

6 

1000 

low 

low 

^h 

55 

0.66 

5.6x103 

0.62 

Alouette I, low 
latitude, night 

7 

1000 

high 

fn> 

fh 

55 

1.49 

i 

2.8xl0‘t 

'1 

1 

; 0.62 

Alouette I , low 
latitude, day 

8 

1 1000 

low 

high 

fh 

5 

0.66 

! 

i 

5.6x103 

i 1.03 

1 

1 

! Alouette I, high 
. latitude, night 


Note: The subscript s denotes values of the quantities appropriate to the location of a hypothetical 

sounder. 













Figure 5.- Model profiles of plasma frequency, f^^, and gyrofrequency , fj^ 
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Lgure 11.- (a) versus h and (b) dh versus f for c£ 

Cumulative delay curves for 40-km laminations are also shown 
part (b) . 
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Figure 12.- Computed ordinary, extraordinary, and Z traces for case 1. 
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Figure 13.- Computed ordinary, extraordinary, and Z traces for case 2. 
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Figure 14.- Computed ordinary, extraordinary, and Z traces for case 3. 
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Figure 15.- Computed ordinary, extraordinary, and Z traces for case 4. 
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Figure 16.- Computed ordinary, extraordinary, and Z traces for case 5. 
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Figure 17.- Computed ordinary, extraordinary, and Z traces for case 6. 
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Figure 18.- Computed ordinary, extraordinary, and Z traces for case 7 
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Figure 19.- Computed ordinary, extraordinary, and Z traces for case 8 
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Figure 20.- Computed ordinary traces for high £„ and high fj, 
profiles for 6 = 5°, 10°, 15°, 30°. 
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